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Abstract We show how leading radiative corrections can be 
implemented in the general description of h —> At decays by 
means of Pseudo Observables (PO). With the inclusion of 
such corrections, the PO description of h —> At decays can be 
matched to next-to-leading-order electroweak calculations 
both within and beyond the Standard Model (SM). In par¬ 
ticular, we demonstrate that with the inclusion of such cor¬ 
rections the complete next-to-leading-order Standard Model 
prediction for the h —> 2e2fi dilepton mass spectrum is re¬ 
covered within 1% accuracy. The impact of radiative correc¬ 
tions for non-standard PO is also briefly discussed. 


1 Introduction 

The decays of the Higgs particle, /z(125), can be charac¬ 
terised by a set of Pseudo Observables (PO) that describes, 
in great generality, possible deviations from the Standard 
Model (SM) in the limit of heavy New Physics (NP) [1], 
The Higgs PO are defined from a momentum expansion 
of the on-shell electroweak Higgs decay amplitudes. More 
precisely, the PO relevant to h —? At are defined by the mo¬ 
mentum expansion around the physical poles (due to the ex¬ 
change of SM electroweak gauge bosons) of the following 
three-point correlation function 

(0\&{j?(x),JZ(y)M0)}\0), ( 1 ) 

where (x) are generic leptonic currents. This expansion 
encodes in full generality the short-distance contributions to 
the decay amplitudes in extensions of the SM with no new 
light states [1]. However, in order to compare this ampli¬ 
tude decomposition with data, also the long-distance contri¬ 
butions due to soft and collinear photon emission (i.e. the 
leading QED radiative corrections) must be taken into ac¬ 
count. 


Soft and collinear photon emission represents a universal 
correction factor [2,3] that can be implemented, by means of 
appropriate convolution functions (or, equivalently, shower¬ 
ing algorithms such as those adopted in PHOTOS [4], PYTH- 
IA [5], or SHERPA [6]) irrespective of the specific short- 
distance structure of the amplitude. 1 In this paper we illus¬ 
trate how this works, in practice, in the h —> At case. 

We focus our analysis to the h —> 2e2fi case, that is par¬ 
ticularly interesting for illustrative purposes: the effect of ra¬ 
diative corrections can be implemented by simple analytic 
formulae, allowing a transparent comparison with numeri¬ 
cal methods. As we will show, the inclusion of the universal 
QED corrections is necessary and sufficient to reach an ac¬ 
curate theoretical description of the Higgs decay spectrum, 
that recovers the best up-to-date SM predictions in absence 
of NP. 


2 QED corrections for the h —)• At dilepton spectrum 


In this section we describe how leading QED radiative cor¬ 
rections affect the dilepton spectrum of h —> At decays as¬ 
suming a generic PO decomposition of the amplitude. As 
anticipated, we focus our discussion to the case of two lep¬ 
ton pairs with different flavor (h —> 2e2fi) and, more pre¬ 
cisely, on the double differential lepton-pair invariant-mass 
distribution 


F(m ee ,mnn) 


d 2 r(h ->• 2e2n) 

dm«,d/n MiU 


( 2 ) 


The emission of soft and collinear photons leads to in¬ 
frared (IR) divergences in the h -4 At spectrum. The full 
structure of such divergences is rather complicated. How¬ 
ever, as we have checked by means of an explicit calculation 


'For a discussion about the implementation of universal QED correc¬ 
tions in a general EFT context see also Ref. [7], 
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at 0(a), such divergences can be factorized in F(m ee ,m ^ lu ) 
and can be analyzed separately for each dilepton system. 
This happens because each fermion current in Eq. (1) car¬ 
ries an overall neutral electric charge. 

Working in the limit of massless leptons, we need to in¬ 
troduce two independent IR regulators for soft and collinear 
divergences. We choose them to be: i) the minimal fraction 
of invariant mass lost by the dilepton invariant-mass sys¬ 
tem; ii) the minimal invariant mass of a single lepton plus 
(collinear) photon (m*). 

We then define the radiator co(x,x*), that represents the 
probability density function (PDF) that a dilepton system re¬ 
tains a fraction sfx of its original invariant mass after brems- 
strahlung for a given x* = 2ml/m^, where is the initial 
dilepton invariant mass (pre bremsstrahlung). By construc¬ 
tion, the kinematical range of x is 

0 < X < Xmax = 1 X* . (3) 

Keeping only the leading terms for (1 — x) « 1 andx* <C 1, 
the radiator is 


®(x,X*) = ft)l(x,X*)0(l — X* — x) -I- 0>2(x,x*)(5(l — x) , (4) 

where 


, . at l 
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n \ 1 — 


log 


2(1 -x)-x« 

X* 


a 
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(5) 


The first term, ft)], describes the real emission of a photon 
such that the lepton pair retains a fraction y/x of its invariant 
mass; the 0-function implements the corresponding IR cut¬ 
off. The second term, ah, describes the events in which the 
soft radiation is below the IR cutoff, as well as the effect of 
virtual corrections. 

We have determined the structure of ft)] by means of an 
explicit 0(a) calculation of the real emission, while Oh has 
been determined by the condition [q dxft)(x,x*) = 1. The lat¬ 
ter condition implies a redefinition of 0(a/n), not enhanced 
by large logs, of the PO characterizing the non-radiative am¬ 
plitude. 

Denoting by men and mo 2 the invariant masses of the 
two dilepton systems before bremsstrahlung, defining fur¬ 
ther x,' = (mj/moi) 2 , it is easy to show that 


d 4 r 

dmoidwo2dx]dx2 


F 0 (mo 1, )«02) (0 (x 1, X1 *) ft) (x 2 , X2* ) , 


(6) 


whereFo(moi,wo 2 ) denotes the non-radiative (tree-level) spec¬ 
trum [1], 

Starting from Eq. (6) we can extract the double differen¬ 
tial spectrum after radiative corrections. To this purpose, we 
first trade xi i2 for m\, 2 , obtaining: 
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From Eq. (7) we then explicitly extract the double differ¬ 
ential decay width by integrating over all the possible physi¬ 
cal moi o 2 combinations, determined by the conditions moi + 
m 02 < )«/, and moi,o 2 > m\.i/yjx mm . In this way we finally 
obtain: 

fm h —m 01 4mim2 


rm r 

F(mi,m 2 )= / __ dm 0 i / 


I xj m^+2m~ 
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xF 0 (moi,m 0 2)ftt ( -^r,xi* ) 

V'«oi J 


co ~^r,x 2 * . (8) 


02 


We stress that the result in Eq. (8) includes both real and 
virtual QED corrections. The latter have been indirectly de¬ 
termined by the normalization condition for ®(x,x*), that 
is the same condition applied in showering algorithms [4], 
As anticipated, this implies a 0(a/n) redefinition of the PO 
compared to their tree-level values (both within and beyond 
the SM). In the context of next-to-leading order (NLO) ef¬ 
fective field theory (EFT) calculations [8,9], this procedure 
provides a well-defined condition for the matching between 
the full EFT calculation of the amplitude and the PO decom¬ 
position. 


3 Comparison with full NLO electroweak corrections 

In this section we present a comparison of the SM predic¬ 
tions for the h —> 2e2fi dilepton invariant mass spectrum ob¬ 
tained using full NLO electroweak corrections [10], and the 
PO decomposition “dressed” with leading QED corrections, 
as described above. 

The complete SM NLO electroweak corrections to h —» 
4£ have been computed in [10], and the results have been im¬ 
plemented in the Monte Carlo event generator Prophecy4f [1 1] 
We have used Prophecy4f version 2.0 to generate 200 mil¬ 
lions weighted events for the recombination mass parameter 
m* = 1 GeV. We have used the default Prophecy4f SM in¬ 
puts except for setting the Higgs boson mass to 125 GeV. 
Prophecy4f adopts the dipole subtraction formalism [12] for 
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Fig. 1 Left: Dilepton invariant mass spectrum in the SM for the It —> 2e2fl decay (full line: PO decomposition “dressed” with QED corrections; 
red and blue bands: complete NLO result from Prophecy4f). Right: Dilepton invariant mass spectrum in the presence of new physics for various 
benchmark scenarios (see text for details). 


the treatment of soft and collinear divergences, and the so- 
called “photon-recombination” is applied. In particular, if 
the invariant mass of a lepton and a photon is smaller than 
m*, the photon momentum is added to the lepton momen¬ 
tum [10] .Asa result,coincides with the collinear cut-off 
introduced in the previous section. 

In Fig. 1 (left) we show the decay distribution as a func¬ 
tion of the dilepton invariant mass normalized to the total 
decay width for h —> 2e2u in the SM (upper plot) and the ra¬ 
tio between NLO and leading-order (LO) predictions (lower 
plot). Shown in solid black is our improved prediction ob¬ 
tained by convoluting the leading order distribution, shown 
in dashed black, with the radiator function as described in the 
previous section. The PO have been fixed to their SM tree- 
level reference values (Kzz = 1, £, = 0 [1]). The Prophecy4f 
predictions within MC uncertainty are shown with red and 
blue bands for and e + e~ invariant mass spectra, re¬ 

spectively. 

We list here a series of conclusions that can be derived 
from this numerical comparison. 


- The spectrum obtained with the PO decomposition of the 
amplitude, “dressed” with leading QED corrections, pro¬ 
vides an excellent approximation (within 1% accuracy) 
to the spectrum obtained with full NLO EW corrections. 

- The effect of the leading QED corrections can be large, 
exceeding 10% in specific regions of the phase space. 
It therefore must be included, in view of a precise data- 
theory comparison, also when fitting beyond-the-SM pa¬ 
rameters. 

- The PO “dressed” spectrum is obtained setting £, = 0 
(i.e. to their LO SM values). The good agreement with 
the complete NLO calculation confirms that the 0{a/n) 
redefinition of the e, is a small effect, with no observable 
consequences for the h —> 2e2fi dilepton invariant mass 
spectrum. 


2 The ~ 2% deviations at the border of the phase space are expected 
due the breakdown of the approximation mu 3> m„ employed in the 
analytic evaluation of the radiation function. 





























4 


4 Implications for New Physics 

As shown in Fig. 1 (left), radiative corrections can be siz¬ 
able and must be included also when going beyond the SM. 
Having demonstrated the validity of our QED improved pre¬ 
dictions to describe such effects, we are in position to ap¬ 
ply the method in the presence of an arbitrary new physics 
contribution to h —> 2e2ji decay as parameterised by generic 
PO [1], As an illustrative example, we consider the impact of 
the leading QED corrections for non-standard values of ic zz , 
£ze L , £ze R . fc'z/j, . and £ ZjUj! . 

To draw some general conclusions we analyse three ben¬ 
chmark points, chosen such that the deviations of the total 
h —> 2e2jj. decay rate from the SM prediction are always 
small, 1 but the impact on the spectrum are quite different. 
The results of the inclusion of QED corrections are shown in 
Fig. 1 (right). As in the left panels, we plot the dilepton in¬ 
variant mass distribution normalized to the total rate (upper 
plot) and the ratio between NLO and LO (lower plot). 

The definition of the benchmarks, and the consequences 
following from the analysis of radiative corrections, are listed 
below. 

- Benchmark I [jc zz = 1.3, = £ Zm , = -0.05 , e ZeR = 

£ zh r = 0.05 (dot-dashed blue)]. 

Here the deviation from the SM point in the Higgs PO 
parameter space is small: this benchmark point is com¬ 
patible with naive power counting in the linear EFT ex¬ 
pansion. As a consequence, small deformations in the 
spectrum are obtained (upper panel) and the relative QED 
corrections are SM-like (lower panel). In this regime, the 
leading QED corrections can be directly extracted from 
the SM result (via an appropriate NLO/LO re-weighting). 

- Benchmark II [ff zz = 0, e Z e L = £'z M/ = 0.26, e ZeR = 
£z/i r = 0 (dotted red)]. 

Here the deviation from the SM point is sizable, beyond 
the naive power counting within a generic EFT (both lin¬ 
ear and non-linear). However, the PO configuration is 
such that the deviations from the SM in the spectrum are 
small. This implies that the relative impact of QED cor¬ 
rections is still SM-like. 

- Benchmark III [rc zz = 0.3, £ Ze , = —0.45 and e z ^ L = 
£ze R = £zh r = 0 (solid green)]. 

In this example we observe a sizable distortion of the 
dilepton shape (upper panel). As a consequence, the rela¬ 
tive impact of the QED corrections is quite different from 
the SM case (a description of radiative corrections by 

3 The dependence of the total rate on the PO can be found in Ref. [13]. 


NLO/LO re-weighting of the SM result would not pro¬ 
vide a good approximation). 

5 Conclusions 

The dominant electroweak corrections to h —> 4/ decays are 
due to the universal soft and collinear photon emission. As 
shown in Fig. 1, these can lead to distortions of the dilep¬ 
ton invariant spectrum of 0(10%) is specific regions of the 
phase space. These effects are of the same order as the ex¬ 
pected modifications from the SM under the assumption of 
underlining linear EFT [ 13], It is then mandatory to properly 
incorporate these corrections in a consistent way both within 
and beyond the SM. 

As we have shown in this paper, this can be achieved 
in general terms within the framework of the Higgs PO. In 
particular we have shown that: (i) the QED corrected predic¬ 
tions for the h —> 2e2ji dilepton invariant mass spectra, with 
PO fixed to their SM LO values, are in agreement with the 
full NLO electroweak SM predictions within 1% accuracy; 
(ii) the QED corrections in the presence of NP can be sizable 
and significantly different from the SM case. 
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